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Abstract 

We address the dynamics of damped collective modes in terms of first and second mo- 

O ' ments. The modes are introduced in a self-consistent fashion with the help of a suitable 

^ I application of linear response theory. Quantum effects in the fluctuations are governed 

by diffusion coefficients D^^. The latter are obtained through a fluctuation dissipation 

, . ■ theorem generalized to allow for a treatment of unstable modes. Numerical evaluations 

'V^ . of the -D^i/ are presented. We discuss briefly how this picture may be used to describe 

5t ! global motion within a locally harmonic approximation. 
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1 Introduction 

The interplay of dissipation and quantum effects has been a long standing challenge to trans- 
port theory, in particular the question of the influence of damping on quantum tunneling. 

*e-niail: hhofniann@.physik.tu-niuenchen.de 
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Convincing answers have only be found in the 80'ies by clarifying demanding problems of 
principal nature with the help of functional integrals, see e.g. ||l|, 0,0, [Q, 0. In these 
approaches one first derives effective actions for global motion from which one may then de- 
duce expressions for the decay rate of a meta-stable state. In contrast to this procedure, we 
like to look at the problem from the point of view of transport equations such as the one of 
Kramers [^ but extended to include quantum effects. By its very nature, such a transport 
equation reflects local properties of collective motion as displayed through the propagation in 
collective phase space. In the following we like to concentrate on situations where damping 
cannot simply be treated perturbatively. For the case of stable modes, suggestions of possible 
quantal equations have been made for instance in |^ and p|. The possibility of generalizing 
to unstable modes has first been described in 0. 

Different to those papers mentioned previously which deal with dissipative tunneling and 
which exploit functional integrals, we like to base our analysis on systems where collective 
variables are introduced self- consistently. For Hamiltonian dynamics, this implies that in the 
Hamiltonian 
@tothamdef 

'H = i/intr + -f^coupl + -f^coU (l) 

for the total system both the "collective" part ifcoii as well as the coupling term are bound to 
be functionals of the dynamical variables of the intrinsic system. One prime example is the 
electron gas for which D. Bohm and D. Pines have been able to deduce such a Hamiltonian 
in |TD|. Another example is nuclear collective motion for which an adapted version of the 



Bohm-Pines method can be found in |11], see also |12|. Although later on we are going to 
concentrate on systems for which such a Hamiltonian ([1|) exists it may be mentioned that our 
method of accounting for self-consistency in the fluctuating force can be generalized to other 
cases, the most prominent example being a situation in which the dynamics is described by 
the Landau- Vlasov equation with (collision term) (see []T3). For stable modes this method 



see 



may in some sense be considered a generalization of the suggestion made by Landau [|1J 

[0)- 

Evidently problems like the motion from a meta-stable minimum across a potential barrier 

are manifestly non-linear in at least one of the collective degrees of freedom. Moreover, there 

are systems for which the intrinsic sub-system and its coupling to the collective part change 

along the dynamical path, thus prohibiting to linearize this coupling in the bath variables. 

In such a situation, as it is given for example in nuclear fission, the technique of functional 

integration is not useful anymore as the internal degrees of freedom can no longer be integrated 

out. It has been argued that it might be possible to handle such systems by a locally harmonic 

approximation (LHA) (see e.g. |]K[, [|r^, IjT^). In this method it is not necessary to restrict the 



dynamics of the internal degrees of freedom to a system of (perturbed) oscillators. Certainly, 
as an essential intermediate step, there will appear the dynamics in linearized version for 
the collective degrees of freedom. As already demonstrated by Kramers, in lucky cases it 
may even suffice to have information on the dynamics at the potential minimum and at the 
barrier, the latter mode being unstable. In this paper we shall address specifically these 
problems and study the influence of dissipation on average motion as well as on the dynamics 



of the second moments (fluctuations). This will be achieved within the framework of linear 
response theory. The latter does not only allow one to treat average motion (in a linearized 
fashion), it opens an ideal playground for deducing the diffusion coefficients which govern 
the dynamics of the second moments. To be able to treat instabilities as well, it will turn 
out necessary to perform suitable analytical continuations. Manipulations of the kind which 
will be needed here have been used before in ||18[, [|l^, but, as mentioned, not in connection 
with a transport equation and, hence, diffusion coefficients. As we shall see the numerical 
computations of these coefficients will very naturally hint to limits of the LHA at smaller 
temperatures. 

2 Brief review of average motion 

In this paper we would like ourselves to restrict to a situation in which the dynamics of interest 
can be described by the one quantity F. For average motion this means looking at the time 
evolution of the expectation value (-F)t- To start with a simple situation let us suppose we 
are given a Hamiltonian of the type 
@twobodham 

H^'^ = Ho + \fF (2) 

with some coupling constant k, without bothering for the moment how it has come about. 
Later on we are going to describe how this form may derived from more general considerations 
and how the k may be deduced from static properties of the system, in which case it will 
depend on the actual state of the intrinsic system. Let us suppose further that we are 
interested in harmonic vibrations. As may be well known, their frequencies and strengths can 
be deduced from the following response function 
@chicoll 

( \ X(^) fo\ 

^-^^(") = iTk^) ^^^ 

obeying the following definition 

@qext 

5{F)^ = -Xcoii(w)/cxt(^) (4) 

The S{F)i^ measures the Fourier transform of the deviation of {F)t from the equilibrium value 
(-F)o, as it follows from the response of the system in linear order to an external field fey±{uj), 
with the coupling being given by 5H = fextF. The x{^) is the Fourier transform of the 
response function defined by 
@defresmic 

X{t -s) = 0(t - .)^tr (pl\F\t), F^{s)]) (5) 



with the Pqs representing the equilibrium density operator associated to the Hamiltonian ^o- 
Notice please that this function is governed entirely by intrinsic properties, different to the 
Xcoii(^) which contains information on collective motion. 



We should expect average dynamics to be closely related to the mean field approximation. 
Indeed, for the latter one readily verifies the following Hamiltonian 
@hammeanf 

^S = Ho + {Q- Qo)F (6) 

where the abbreviation 
@selfcons 

k{F)t = Q-Qo (7) 

has been used. One may now simply apply the Clausius-Mosotti construction to derive the 
form (|) from (^) and (0). 

For the following it will show convenient to add to the Hamiltonian H^l of (|^) a c-number 
term and thus to introduce a H{xi,pi, Q) like 
@hamexp 

H{x.,p., Q) = H{xup„ Qo) + (Q - Qo)Hs^^,P^, Qo) + liQ- Qof ( ^'^^y'^^ Qo • (8) 

Identifying H{xi,pi, Qq) = Hq and the expectation value on the very right with 
@def coup con 

, -1 _ / d^H{xi,Pi,Q) ^ 
~^ "^ dQ' ^^° ^^ 

one sees that the total energy of the system i^tot = (-^^mf ) ~ 2 (-^) (-^) ^^^ simply be expressed 
as the expectation value {H{xi,pi,Q)) of the one-body Hamiltonian H{xi,pi,Q). 

This discussion has demonstrated the intimate relation between a Hamiltonian of the form 
(H) and the one of (^) involving a schematic two-body interaction. On the level of the mean 
field approximation the latter becomes equivalent to the former. However, the H^'^^ of (0) 
allows one to treat fiuctuations in F, or what will turn out more convenient in Q — Qq, about 
its average value. Indeed, later on we will be exploiting exactly this feature when we are 
going to quantize the "collective variable" Q. Conversely, as often mean field approximations 
are more easily accessible numerically than solutions of the full Schrodinger or Heisenberg 
equations (of the many body problem), one may start from such a H{xi,pi, Q) and construct 
the effective two body interaction appearing in H^"^' . 

A prominent example from nuclear physics is given when one tries to simulate the mean 
field approximation via the deformed shell model. In that case the Q is introduced as a 
c-number variable specifying the shape of the nucleus. The latter then appears both in the 
single particle potential as well as in the liquid drop energy. To account for the latter is 
necessary because otherwise one would not be able to calculate reliably the total energy of 
the system. This can be done with the help of the Strutinsky renormalization. How one 
may then proceed to construct a Hamiltonian of the type H{xi,pi,Q) has been worked out 
in detail in ||20[| . 

It is more than conceivable that one might start from even more microscopic theories like 
the one of Hartree-Fock involving effective forces or the Relativistic Mean Field theory. An 
equation like (0) would then appear as a constraint for possible variations in the deformation 



of the self-consistent field. After appropriate linearizations and after applying time dependent 
perturbation theory in suitable manor one would end up in a form like (|^), provided that one 
is interested in the change of the expectation value of the same operator F which specifies 
the variation in the density. In general the response "function" would have tensor character, 
of course. 

In this context one may mention examples worked out for Landau Fermi liquid theory 
in [21 1, |13|. Looking at general modes of the single particle density in phase space for the 
associated response functions Xcoii(^) forms similar to the one given in (||) have been deduced. 
There it might only happen that the " coupling constant" k may depend on the wave vectors 
of the various modes. 

2.1 Extension to finite excitations and large scale motion 

In the following we would like to concentrate on the example for which a Hamiltonian 
H{xi,Pi,Q) is given, which depends parametrically on a c-number Q{t) representing the col- 
lective degree of freedom. It is this example which is developed farthest theoretically and 
which still can be handled numerically even for such complex situations as nuclear fission 
[p^ , [p3[] . As described in various papers (see e.g. [^ and [|l^], in particular), with such a 



Hamiltonian it is possible to derive an intrinsically consistent theory for collective motion of 
large scale provided there is a clear separation of time scales. 

Let Tcou be a relevant measure for a collective time scale and r the one for the remaining 
degrees of freedom, henceforth called the intrinsic ones. For time lapses St with t < 6t '^ 
Tcou one may describe collective motion within a (locally) harmonic approximation. This 
is achieved by exploiting the expansion (||) and identifying simultaneously the unperturbed 
density operator p^g as the one corresponding to the Hamiltonian H{xi,pi,Qo). This density 
operator may then either be a function of temperature or entropy depending which distribution 
we want to use for parameterizing the equilibrium. For the sake of simplicity, in this paper 
we will choose the canonical one later on, but we may refer to |T^ for a discussion of the more 
general case. Thus we may write for the density Pqg = Pqs{Qo, Tq), with the Qq and the Tq to 
represent those values which the "macroscopic" variables Q,T have during the time interval 
St. 

As shown in [53] and []12| all the steps discussed above for the genuinely harmonic case 



may then be taken over. We must expect, of course, that all quantities will depend on the 
actual quasi-static state of the system which is being specified by the pair of variables Qo, ^o- 
This is true in particular for the coupling constant k which can be written in the form 
Icoupcon 

+ x(0) = C(0) + x(0). (10) 



_^-i= d^E{Q,So 



Here E{Q., Sq) is meant to represent the internal energy of the system, with the entropy 5*0 
being understood to be calculated for given Qo,Tf). (It is possible, of course, to rewrite this 
formula in terms of derivatives of the free energy). The x(0) stands for the static response and 
the C(0) is introduced as a short hand notation for the local static stiffness. To derive this 



relation, together with the form (||) for the collective response function, one needs to apply 
similar arguments as mentioned before but has to care for possible changes in the quasi-static 
properties of the system. In doing so one may be guided by the fact that within the harmonic 
approximation (for the motion in Q) entropy can be assumed to be constant (to leading 
order), as it will change with velocity only quadratically Using this hypothesis, in formal 
sense the derivation of the form of the collective response XcoiK"-^) follows identically the one 
given above and which ended in (^, provided the adiabatic susceptibility x'^'^ is identical to 
the static response, x^*^ = x(0); for a general discussion of the relevance of this property see 



2| and |12 



We should like to note that for the general situation away from the potential minimum, the 
Qo in (1^ is to be replaced by the Qm defining the center of that oscillator by which the "true" 
static energy E{Q,Sq) is approximated in the neighborhood of the Qo of (H) (see e.g. |jl2[). 
Furthermore, we want to mention that the discussion to come shall be restricted to situations 
of negative coupling constants, which for the case of an instability will imply the | C(0) | to 
be smaller than x(0), see below. Looking back to the two body interaction introduced in (0) 
this is to say that we like to restrict ourselves to effective forces being attractive. For the 
nuclear case this is known to correspond to motion with neutrons and protons in phase. 

2.2 Reduction to single modes 

The poles of Xcoii(^) define possible excitations of the system, more precisely those which 

involve dynamics in the quantity (F) and which are caused by that perturbation SH which 

is proportional to the F itself. Considering all of them would in the end imply strict non- 

Markovian behavior in {F)t- Eventually the strength distribution is dominated by only one 

mode, in which case the Xcoii(i^) may be replaced by the response function of the damped 

oscillator 

@resosc 

1 1 

xLH = - .jF 2 I ■ F 2 (11) 

M^ uj'' + zujL — w^ 

with M^Y = 7^ and M^vj"^ = C^, where 7^, M^ and C^ stand for the coefficients of friction, 
inertia and stiffness for the (local) motion in S{F)t, a feature which is easily verified from (Q). 
For several reasons it may be advantageous to introduce the coordinate q{t) = Q{t) — Qm for 
which the local equation of motion writes 
@eomaverq 

Mq{t) + 7g(t) + Cq{t) = -q,^t{t) (12) 

This transformation is readily performed considering (|^), with the Qo being replaced by Qm, 
of course. It implies to work with the gg-response function 

@qqresp 

5q{uj) 



k Xcoll(^) = Xqqi.^) 



(^(-goxt(^) 



(13) 
gcxt=o 



instead of the Xcoii(i^), with a similar relation leading from the X^cl'^) ^o ^^e oscillator response 
Xosc(i^) for q-motion. The latter is then specified by the transport coefficients 7 = A;~^7^, M = 



k-'^M^ and C = k-^C^. 

Whether or not such an ideal situation is actually met depends largely on the value of 
the coupling constant k. In fair approximation one may say it to be given whenever in (|T0|) 
the static stiffness C(0) is small compared to the static response x(0)- In such a case the 
strength of all the solutions of the secular equation 1/k + x{^) = will concentrate in the 
low frequency domain. Therefore it is important to realize that the static energy may depend 
strongly on temperature. For nuclear physics, for instance, the C(0) drops dramatically above 
ksT = 2 MeV where the influence of shell structure has diminished (see e.g.|p5|). 

In the more general case the strength distribution of the collective modes, which is deter- 
mined by the imaginary part Xqq{^) of the gg-response function, will spread over several (or 
even many) modes. In this case the equation of motion for q{t), which may be obtained from 
@qeomnm 

{Xqq{^)y^ Qi^) = -gcxt(^) (14) 

after a Fourier transform back to time, will contain terms being non-local in time which 
cannot be expressed simply by a sum of derivatives of q(t) up to second order. If one still 
wants to stick to differential form one needs to apply a further approximation. For instance, 
one may simply expand the {xqq{^))~ in (HH) to second order in u [^]. Another possibility is 
found in "reducing" the Xqqi^) to the Xosc(^)- By this we mean to replace the former by the 
latter in a certain range of frequencies. In practice this may be done by fitting the imaginary 
part of Xosc('^) to that peak in Xqqi^) the corresponding mode of which one wants to treat. 
Notice please that both of these approximations require that only frequencies are relevant 
which lie in some limited regime. In this sense the second one may be considered superior to 
the former. In the past, various computations of transport coefficients have been performed 
which lead to values being in fairly good agreement with experimental experience, see [^, 
2| 



2.3 Unstable modes 

The oscillator response is analytic in the stiffness C. As we will see below, the latter becomes 
negative whenever the static stiffness C(0) becomes negative, which happens in the region of 
potential barriers where the system is unstable, in the sense that one of the two fundamental 
solutions for average motion begins to grow exponentially. Such a situation can still be 
handled within the locally harmonic approximation as long as the forces do not change too 
rapidly with the coordinate. Recall please that at some given Qq the harmonic solution is 
meant to portray the true one only within a limited time interval 6t anyhow. This growth of 
the time dependent solution reflects itself in the feature that for negative C one of the poles 
of the collective response moves into the upper part of the frequency plane. (Notice please 
that the instability we have in mind still leaves both inertia as well as friction positive). This 
does not cause any problems when dealing with the functional form of the oscillator response 
as given by (pUD, but it may well be so for more general expressions. For instance, one may 
want to write the response function in terms of integral representations, like 



©rftodrfst 

Xf,u{^) = — 77 forlmcj>0. (15) 

Jc vr \l — uj 

In this case some caution is required in defining the contour C. For stable systems it is taken 
along the real axes with the poles lying below. This choice is required by causality. If the same 
principle is still to be valid for the case of an instability the contour C must still be chosen to 
be above the one singular pole mentioned before. Actually, besides this "retarded" function it 
may become desirable to have an " advanced" response function as well, for which the contour 
then has to be modified to lie below the poles which correspond to motion backward in time. 
More details of these formal techniques are described in Appendix A. 

2.4 Implications from sum rules 



The integral representation (|T^) may serve as a starting point for deriving and exploiting sum 
rules, even in the case of instabilities. Generalizing the common definition to 
@sumruldef 

Sin) ^ f ^^n.^^^ ^ (ig) 



^ Jc TT 

the same relations are recovered for the unstable case as those known for the stable one, if 
one only observes the appropriate choice of the contour. For example, for n = —1 the value 
of the integral in (|T6|) gives the static response, Sj~^^ = Xnv{'^ = 0), while for n = 1 one may 
employ the usual arguments to obtain 
@suruone 

Ic n ^"^ ' ^^°o ic vr 1- 



^ = / v"C(f') = Jto -i v™ = - Jl"i-^-(-) ■ (") 



One way of getting the equality in the middle is to involve the geometrical series for the 
1/(1 — (f2 — ie)/uj)^ which is possible whenever the 5*^"^ remain finite for n > 0. The latter 
restriction is not really necessary. One may simply replace the Vt in the first integral by 
VL = uj{l + Q/uj) = [1 — (Q/uj) — {Q/uj)^/{l + Q/uj)]^^ and observe that for u; — i> cxo the last 
term in the brackets gets less and less important. Then one only needs to have the SJ^J vanish, 
which in our case is no problem as later on we are going to apply ( p!7| ) only for the diagonal 
case. Please observe that for the oscillator response (pH]) the validity of the expression on the 
very right follows directly from inspection. 

The sum rules for ra = =pl are intimately related to two of the three transport coefficients. 
This can be inferred from the values these two sums take on for the oscillator response, namely 
@sumrulval 

SL'^ = ^ Sill = l^ with 5^ = 0. (18) 

For later purpose we may note here that these results remain unchanged when the oscillator 
is modified by introducing a frequency dependent friction 7 = 7(0;). 

As the oscillator response is derived from the original, microscopic expression given by 
(^ together with (|T3]), the sum rules may be used to derive inequalities for the transport 
coefficients. How these inequalities come about may best be seen writing the full response 
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function in terms of a pole expansion like 
©resoscsum 



'^»<'^' = -?2iio:to-:^ (I'" 



a 



with C(j^ = ±£^Q — i^a/'^- Then one may deduce results like those given in ([T8| ) for each term 
of the sum over a. It is then evident that the value of the 5*^^^^ for each term cannot be larger 
than the corresponding value of the total sum. 

Let us look at the inertia first. The form (|^) implies lim^^oo'^^(Xcoii(^) — x(^)) = 0- 
Therefore we get 
@sumqqin 

k-'sfi= r-x'Lmn= r-x"mn = U[F,F]) (20) 



oo 



/vcoiiv / / /v v--/ -- f- 

vr J-oo vr n 



and consequently 
@sumrulinert 

^<5«^f(|F.F,).-i- (21, 

which means that the inertia for the one mode we have chosen to pick cannot be smaller than 
the mo given by the microscopic value of the energy weighted sum. In nuclear cases this mo 
can be seen to correspond to the inertia one would get in the liquid drop model for irrotational 
flow (see e.g. |20[). 

Next let us look at the stiffness. Here the situation is slightly different as the mode we 
like to look at explicitly may become unstable. However, assuming that all the other ones 
remain stable one may simply repeat the arguments from before such that for S^~^^ = X(?g(0) 
one gets 
@ineqstiff 

c - "'"'^'^^ - 1 + X(0) " C(0) x(o) + c(o)- ^''^ 

These expressions for Xgg(O) are obtained after observing (|]), ( 1l3D and (p^. The last equation 
shows the X(?g(0) to have the same sign as the static response C(0), from which observation 
it follows that C will become negative as soon as C(0) is negative. In these arguments we 
have assumed the static intrinsic response x(0) to be positive, as there is no reason to expect 
the intrinsic system to be unstable. As mentioned before, we want to stick to the case where 
this value is larger than the absolute value of the static stiffness. The inequality allows one 
to prove that effective stiffness cannot be smaller than the static one: 
@ineqstiffres 

C > C(0) (23) 

For stable modes this is immediately clear, for unstable ones it is useful to make a detour and 
introduce the absolute values C = — \ C \ and C(0) = — | C(0) |, in which way one gets from 

m 

@ st iffb oundf einv 

111 1 

^ I ^/ X I 1 ^, , I , > , ^, . , . (24) 



C I - I C(0) I 1- I C{0) I /x(0) - I C{0) 



The inequality (^) expresses in clear fashion that in general the effective stiffness C will be 
different from the static one, C(0). The reason for this difference is found in dynamical eSects 
contributing to the effective potential. It is interesting to note that at an instability these 
effects actually lower the absolute value of the stijfness, opposite to the case in stable modes. 
In this sense the dynamics weakens the instability by c?e-creasing the conservative force. 

3 Inclusion of dynamical fluctuations 

In the more general sense our treatment so far has not gone beyond a mean field approxi- 
mation. It is true, however, that the dynamics in the basic quantities {F)t or q{t) may be 
damped. Microscopically this damping implies that we need at least to consider effects which 
are not contained in a pure Hartree (or Hartree-Fock) approximation. However, it still suffices 



to stick to a description in terms of self-energies in the single particle degrees of freedom ||T2 
and the q{t) still simply is a c-number. Evidently the latter restriction has to be given up 
whenever one is interested to describe processes where fluctuations in the collective degrees 
of freedom become relevant. Take the example mentioned before, the decay out of a potential 
minimum. One knows from the pioneering work of Kramers how in the realm of classical 
physics the dynamics of fluctuations across the barrier can be treated and how they influ- 
ence the decay rate. In E^ it was shown how the quantal extension of this decay rate can 



be obtained within a description based on transport equations. Some of the ideas behind 



this derivation were borrowed from the method developed in [0. Unfortunately, also this 
method was again based on functional integrals for which reason it can be applied only to the 
schematic cases mentioned in the introduction. In this section we thus would like to proceed 
describing the basic steps needed for a formulation in terms of transport equations, with an 
elaborate discussion to come in the last two sections of the paper about the essential ingre- 
diences of this transport equation. By those we mean these quantities which are intimately 
associated to fluctuations, namely the coefficients which parameterize the diffusion process in 
collective phase space. 

3.1 Quantization of the collective variables 

As the first step we need to render the collective degree of freedom to become a fluctuating 
variable. Often people are tempted to do this simply by adding to the Hamiltonian (^ a 
kinetic energy term like 11^/ (2m), with some unknown mass parameter m to represent the 
collective inertia and 11 being the momentum conjugate to the coordinate, which then is to 
be considered an operator Q. It is not difficult to convince oneself that in this way self- 
consistency would badly be lost; for a detailed discussion see ||12|. This feature cannot be 
cured by playing with the potential part, which means to say by adding some function of 
Q — Qq. Rather one needs to introduce a momentum dependent coupling. 

Starting from the effective Hamiltonian (|^) in the original particle space, this method 
allows one to deduce the following Hamiltonian for the total system (see ||11||): 
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@totham 

n = H{x,,p,, Qo) + kUF-^{Q- Qo)F + ^ + ^ ^^^, \ Q - Qof (25) 

Comparing with (|l]) one easily identifies the parts for bare intrinsic and collective motion, as 
well as the coupling between both, which here consists of the two terms being linear in 11 
and Q — Qq. In (|25D the F = i [H{xi,pi,Qo),F], the 11 stands for the canonical momentum 
and tuq is the inertia which by (pT]) was introduced to represent the energy weighted sum, 
and which serves here to define the unperturbed collective kinetic energy. The parameter (3 
has to be specified only in connection to the collective momentum, a point to which we will 
come to below. This f3 drops out of any equation of motion which only involves the collective 
coordinate Q. This feature is very similar to the properties we found before for average 
dynamics (as expressed in (F) or in q) which could be traced back to the response function 
(^. It should be no surprise that the same function appears again when the dynamics for 
average motion in q{t) is re-derived from Ehrenfest's equations to (|25D . We may note in 
passing that in the original version of the Bohm-Pines procedure, as it had been developed 
for the electron gas in |l^ , no such additional parameter (3 was needed. The difference to our 
case lies in the fact that we want to apply the method to an attractive two body interaction 
in contrast to the repulsive one of the electron gas (c.f. |TT|). 

The description in terms of the Hamiltonian (^) is not yet complete. As it stands the 
Ti would have too many degrees of freedom and in this sense could not be equivalent to the 
original two body Hamiltonian H^'^^ of (@). Indeed, the Bohm-Pines procedure necessarily 
involves a subsidiary condition which for the present case reads 
@ sub con 

kF = Q-Qo (26) 

It is nothing else but the "operator version" of the self-consistency condition (|^). This property 
and the very fact that it is inherently incorporated into the form (p5|) of the Hamiltonian Ti 
makes it very plausible that the description of collective motion on average is the same as 
before. The new feature is seen in the fact that now we are in a position to treat fluctuations 
in the q in the way similar to what the Hamiltonian H^"^^ of (^ allows one to do for the 
microscopic quantity F. This feature goes along with the fact that in (|25|) no two body 
interaction appears anymore. Rather, its effects are hidden in the terms involving collective 
coordinate and momentum. 

With the Hamiltonian (|25|) one may proceed to derive effective equations of motion for 
the propagation in collective phase space. In this spirit the Nakajima-Zwanzig projection 
technique has been applied in to get the transport equation for the Wigner function of the 
collective density. However, to account for self-consistency not only on the level of average 
dynamics, or on that of the mean field, a non-perturbative method had to be developed for 
treating the integral kernel. Again this was possible within the locally harmonic approxi- 
mation for which it suffices to know the propagators for the consecutive time steps 6t with 
St ^ Tcoii- As for 6t -^ this propagator starts at a definite point in phase space, it was 
argued that a Gaussian approximation would do. Thus one only needs to specify the time 
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evolution of the first and second moments. For this non-perturbative Nakajima-Zwanzig equa- 
tion it was possible to prove that it is in accord with the (quantal) fluctuation dissipation 
theorem (FDT) 0], which in full glory requires to have non-Markovian equations of motion. 
This feature makes it possible to take advantage of the FDT for establishing the approximate 
equations of motion, which for the first and second moments may be taken to be of differential 
form. 

3.2 Harmonic approximation for first and second moments 

For average dynamics we have seen the differential equation of motion to come up through the 
reduction of the "full" response function (0) (or the one of (0)) to that of the oscillator given 
in (|TT|). Evidently, one would like to have the equations of motion for the second moments 
to be of similar type. Clearly, in addition to the coordinate they will involve the collective 
velocity or the corresponding momentum as well. Let us first introduce the latter on the level 
of first moments, i.e. for qdt) = {q)t and Pc{t) = {p)t- It is clear that for the harmonic case 
we are looking at, one expects the following set: 
@diffeqqk 

M-ge(t) = Pc{t) -pc{t) = -Ij^Qcit) - Cq,{t) , (27) 

As the first equation tells one, the momentum appearing here is the kinetic one, called p 
henceforth. It is here where one may benefit from having the parameter (3 appearing in (|25|). 
As shown in [|11| and [0, the choice 13 +k = k'^C allows one to replace the canonical momentum 



n by the kinetic one not only for average motion but also for the fluctuations around it. As 
demonstrated in ^ these fluctuations satisfy the following linear set of equations 
@smsigqq 

J^^M - |:S,,(t) = (28) 



ysmsigqp 



ysmsigpp 



^S,,(t) - ^S,,(t) + CS,,(t) + ^S,,(t) = D,, (29) 



-Spp(t) + 2CS,p(t) + ^Spp(t) = 2Dpp (30) 



—T,pp(t) + 2CT,qp[t) + — . 

The second moments are defined as Sgg(t) = {(t)t — qlif)i and similarly for the combinations 
{q-iP} and {p, p}. It is seen that the structure of the homogeneous parts follow from that 
of (0), for which reason the q^if) etc. satisfy the set ( P^D to ( PD| ) with the inhomogeneities 
Dqp, Dpp put equal to zero. For any finite value of these diffusion coefficients Dgp and Dpp the 
fluctuations will become finite in the course of time. This will be so even for zero initial values, 
as it is the case when the solutions are to represent the second moments of the propagators 
(see above and below). So far the diffusion coefficients have not been specified. It is here 
where we may benefit from the FDT. As is easily verified, these diffusion coefficients can be 
expressed by the stationary solutions of (pSD to (0). Let us first take the conventional case 
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of a positive C for which these stationary solutions are those of equihbrium such that we may 

write: 

Odiffcoef 



^PP— nj^PP Dqp — C^qq M^PP (^l) 



7 1 

M PP ^qp-^^qq ^■ 

In the next chapter we will exploit the (quantal) FDT to evaluate these expressions. After- 
wards we will turn to the case of an instability with C < which will be handled by analytic 
continuation. 

Often one encounters cases of strongly over-damped motion. It is not difficult to convince 
one-selves (see ||T2[ and c.f. 0]) that in this case the set of equations reduces to (|28|) to (|30|) 
@smsigqqovd 

J^^M + 2^^qqit) = 20-" (32) 

with the diffusion coefficient D°l'^ = CE^y-f 

It may be helpful to the reader to add some further comments on the seemingly Markovian 
nature of our basic equations of motion. There can be no doubt that approximations of this 
kind are necessary in one way or other. However, for average dynamics this step to differential 
form does involve less stringent conditions than one would expect on usual grounds. This 
feature is intimately related to the construction of the oscillator response or the definition 
of the transport coefficients M, 7 and C. Look once more at (0) and suppose that the 
microscopic response x(^) can be represented by a Lorentzian, or its dissipative part, rather, 
which often may portray a realistic situation quite well. For such a case the form (0) follows 
without any additional approximation. On the other hand, it is almost evident that the set 
of equations (^81) to (0) for the second moments can hardly be correct for all times t and 
in particular not for very small t. This set supposes that we use it for large times, a feature 
already suggested by our way of calculating the diffusion coefficients: They are chosen in such 
a way as to warrant that the set ( p8D to ( pOD describes correctly the relaxation to the proper 
equilibrium. 

3.3 Treatment of large scale motion 



The information contained in (|27D to (0) may now be used to construct the time evolution 
in global sense. Two possibilities come to one's mind. First of all and as indicated before, one 
may take the first and second moments to define a Gaussian. The latter may be constructed 
such as to represent motion within the time lapse 6t = t2 — ti and which at t = ti starts from 
some point Qi, Pi in phase space. For a sufficiently short interval 6t the widths S^i.(t) may 
be assumed to remain small such that this motion remains restricted to a narrow region in 
phase space, thus justifying the linearization procedure we started off. As this propagation 
can be performed for any point Qi, Pi, one may in the end sum up all possibihties needed 
to construct the time evolution of any initial distribution and even over time intervals of 
macroscopically large size. That this "propagator method" may be used even to describe 
such complex problems as the decay out of a potential minimum has first been demonstrated 



m 
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Numerically this method turns out to be time consuming as the "summation over all 
possibilities" requires integration over the phase space not only for the first step but for each 
following intermediate step as well. It so turns out that it is simpler to translate it to the 
picture of a Langevin force. There one follows individual trajectories and constructs the 
final distribution by repeating this latter procedure for a suitably large manifold of points 
representing the initial distribution (in global sense). These trajectories are constructed in 
such a way that for the elementary time step 5t = t2 — ti one first moves along the average 
trajectory, say by solving the set (^Tf) , to find the point representing the actual position at 
time ^2 by Monte Carlo techniques. This sampling is to be performed by taking into account 
the size the fluctuation attains at t = ^2 according to the fluctuating force, which in turn may 
be determined from the information contained in the diffusion coefficient. 

4 Diffusion coefficients for stable modes 



To evaluate the diffusion coefficients from (|3l| ) we need to calculate the values the fluctuations 
would take on if the (collective) system were in equilibrium locally. The fluctuation dissipation 
theorem implies 
@eqaverqq 

duj f huj" 

— coth — - 
27r \2T 



^'-'=^/^^°^hU^P^'^MH- (33) 



Here the indices i/, /i are supposed to represent q and p. Exploiting the relations Xqp{^) = 
iMuJXqq{uj) and Xppi^) = {M^YXqq{^)y ^^^ fluctuatlous can be expressed by the one response 
function Xqq- The different behavior of q and p under time reversal causes the cross term to 
vanish, i.e. E^^ = 0. Formally this is reflected by Xqp being even in uj. It is for this reason, 



that no diffusion coefficient appears in (|28|). 

Actually, the values given by (P3| ) are not yet suitable for our purpose. It was mentioned 
before that the equations (^) to ( pOf ) should be consistent with those for average motion 
given in (pT]). As the latter are derived from the oscillator response we are forced to take the 
latter also when evaluating the fluctuations from an expression like (P3|). In other words, the 
diffusion coefficients ought to be calculated from relations like (|3TD but where the E^*^ are 
replaced by those for the oscillator E°'JJ^, 
@dkinqqpp 

(34) 



(35) 





D,,= 


' yiosc 

M PP 


and L 


v = 


CE°- - 


1 

M 


yiosc 
^pp } 


and the equilibrium fluctuations are 


flnally to be 


calculated from 




@eqaverqqosc 




^qq - 


^ f duj 
h — coth 
J 2n 




x'LM 






and 
















@eqaverpposc 

















Eo;c.M^./^coth(^).^., 



(36) 
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In doing so we keep the E^^ = untouched, as it results from a fundamental symmetry. 

4.1 Limiting cases 

Before turning to the general case let us look at the limits of high temperature and of small 
damping. The former is established if the argument of the hyperbolic cotangent is small, 
i.e. huj/{2T) << 1, for all frequencies which contribute to the integrals. Quite generally, 
one finds: S^^ = Txufii^ = 0) which for the oscillator leads to the classical equipartition 
theorem: 
@equipart 



2M 2 2 

As a consequence, for the diffusion coefficient one obtains the classic Einstein relation 
@einstein 

Dpp = -fT with Dqp = (38) 

Notice that for this result no assumption was necessary about the strength of the friction 
coefficient, besides the fact that 7 will influence the range of frequencies which must be 
considered small compared to T. 

In the limit of zero damping one flnds a similar result, only that the temperature T has 
to be replaced by an effective one, namely 

r-(^) = ^cothg) (39) 

with w^ = \C\/M. This is easily verifled letting the friction 7 in the oscillator response (pT|) 
go to zero. Then the dissipative part of the response function can be written as 
@disoscsd 

hm yL(^) = 7^ (Siu ~w)- 6{uj + U7)) (40) 

The resulting version of the generalized Einstein relation, i.e. Dpp = 7T*, has been used in 
[^ . It allows for a nice and simple evaluation of the correction terms to the classical relation, 
namely Dpp ^ 7^(1 + (l/3)(tJ7/2)^(l/T)^. However, as we will demonstrate below, the form 
Dpp = 7T* approximates the correct formula only for extremely small friction. 

4.2 General case 

To evaluate the integral in ( P5| ) one commonly expands the hyperbolic cotangent into the 
uniformly convergent series 



Dseriecoth 



, fhiu\ 2T "^ 1 ^ 2ttT , , 

coth — = — y = —— 41 



-00 



and applies a general relation (see ( |A.5|) ) which allows one to evaluate the response function 
through its spectral density, the dissipative part; such expressions are briefly discussed in 
Appendix A. For the terms with positive n, the integral leads to the retarded response function 
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calculated at the corresponding pole of the coth, namely 

@contcaus 

.-^^,-X%mB) (42) 

c TX uo — i\n\'a ^ 

For negative n, one may use a relation between retarded and advanced response to find 
@contacaus 



d^ xL(^) 



{xlMn\Q)) (43) 



/c TT a; + z|n|0 

Notice please that in the classical limit it is only the term with n = which survives and 
which leads to Txufi{^ = 0). The quantum corrections are caused by those contributions to 
( ^TD having a finite n. 

A straightforward calculation then shows the diagonal case elements S^^ to be given by 
@eqqqserie 

^7, = T [x,,{^ = °) + 2 E X,M^)) (44) 

In this expression we have used the fact that the response function x^l^l is real along the 
imaginary axis. 

To evaluate the diffusion coefficients we need to apply (|^) to the oscillator response. 
Unfortunately, the resulting series diverges logarithmically for the momentum fiuctuations. 
Actually, such a behavior can better be read off from the integral (^) after recognizing that 
the dissipative part of (0) decreases with uj only like uj~^. To cure this problem one needs to 
regularize the integral. In []^ the latter was simply cut off at a certain upper limit. Here we 
would like to apply the somewhat more elegant method of the so called Drude regularization 
(see e.g. 0). It implies replacing in (^) the friction coefficient by a function of the form 
@drudegom 

Properties of the resulting modified oscillator response function are discussed in Appendix 
|B[ In ( ^51) the additional parameter wd appears, the appropriate choice of which depends on 



general properties of the system. In the present work we shall assume the time scale l/tzr^, 
to be well separated from that of collective motion and will examine how in a physically 
reasonable regime this parameter infiuences the diffusion coefficients. 

For the S°!f one finds the following expression: 
@eqppserie 

J:;i'izuD) = TM h + 2 f; [C + nO^itnQ)] Xosc(m0) j , (46) 

A common procedure to evaluate further both (^) and ( PB| ) is to make use of the digamma 
function "^{z) = -^InT^z). It allows one to sum up series over products of pole terms, for 



instance in the form (see e.g. [^ p^] ) 
@difftwopsi 



^iin + y)in + z) y-z 
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The oscillator response with the constant 7 replaced by the one of ( ^5]) has three poles u^'^' 
(see Appendix ^. Thus one may rewrite it as a sum of three terms where each one factorizes 
like the left hand side of (^T]). After some manipulations one then ends up with 



ysigqqsym 



CE»(.„) = r (1 + ^ g (.,„_j:;^(t;r- .....) * 1 1 + '^ 1 1 («) 




and 
©sigppsym 






where the convention lj*^-'^'^^ = cu^^'^ is used. In (48) also the gg-fiuctuation has been evaluated 



with the Drude regularization, mainly to allow for a clear comparison to the other case where 
this regularization is really necessary. It is clear from (^51) that for large wd the result (^) 
will turn into the one calculated for constant friction. How rapidly this happens depends 
on the values of the parameters involved. This is demonstrated in Fig.l a by numerical 
calculation. It is seen that large values of wd are needed for small values of temperature and 
large damping before the asymptotic limit is really reached. It may be seen as well that even 
for comparatively small zud the difference to this asymptotic limit is quite small. (It should 
be clear from above that wd must not be chosen smaller than the typical frequencies of the 
original oscillator.) Here and in the following it turns out convenient to measure frequency in 



units oi tu = J\C\/M and the degree of damping by the dimension- less ratio rj = 7/(2Mg7). 

Different to the gg- fluctuations, those of the momentum would diverge for large zud- 
Therefore, one must expect them to exhibit a stronger dependence on the Drude frequency 
vud than seen in the gg— fluctuations. This is demonstrated in Fig.l.b. Indeed, for large 
7] this dependence is quite sensible. We may notice that the diffusion coefficient increases 
with wd, although its classical limit is independent of wd. In this sense it can be said the 
quantum corrections to increase with vud, too. The standard values chosen for the figures 
presented here, namely zud ~ 10 u;, are probably very reasonable. Fortunately, in this region 
the dependence of the diffusion coefficient on vjd is not very strong, the value of Dpp changing 
by 30% if wd is increased by a factor 2. 

Let us turn now to the dependence of the diffusion coefficients on temperature T and on 
the effective damping rate t]. Although the transport coefficients for average motion, M, 7 
and C, must be expected to vary with these parameters as well, this feature will be neglected 
here. All computations will be done in Drude regularization. If not stated otherwise, the 
value of zud will be chosen such that the ratio ( = zu/wd equals 0.1. To simplify notation 
we want to use the same symbol as if we would use the mere oscillator, which is to say that 
in expressions like the ones on the left hand side of (|48D and (^) the argument wd will be 
left out. We shall also omit the index osc at times. 

Let us begin looking at the equilibrium fluctuations (see also 0), flrst in the coordinate 
g. Rather than showing the E^^*^ itself, we plot CE"^^, normalized to its value as given by the 
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quantal ground state (Fig.2.a). Fig. 2. a shows the temperature dependence for different values 
of 7]. One finds the quantum effects to disappear both with increasing T and rj. The first effect 
is expected, although it may be somewhat of a surprise how quickly the high temperature 
limit is reached. The second effect causes the E^^ to reach this limit the faster the larger the 
damping. While behaving similarly with increasing T, the fluctuations in momentum differ 
in their dependence on rj: They increase with increasing damping. This is demonstrated 
in Fig.2.b. Here T,p^p/0.5Mzu (or equivalently Dpp/O.bjw) is plotted versus temperature for 
different r/. This figure clearly shows how much the quantum effects are enhanced by friction. 
Take the regime of T < tu where the change with T is quite small. For large rj the fluctuations 
exceed several times the values they have for undamped motion. 

Next we study the diffusion coefficients. They are plotted in Fig. 3 by the dashed lines 
as function of rj for different temperatures, the latter being normalized to w. In Fig. 3. a we 
show the diagonal one, divided by its high temperature limit, Dpp/'jT; Fig.3.b represents the 
non-diagonal one, divided by temperature. Both figures demonstrate that for the diffusion 
coefficients quantum effects increase with damping, and that they may become quite large. 



Take the non-diagonal element plotted in Fig.3.b and recall its definition according to (|3^) . 
It measures the difference between CS^g^^ and j^T^p^p, which vanishes in the high temperature 
limit. But Fig.3.b demonstrates that for small T this difference, and thus Dqp, reaches values 
which are several times this unit. Moreover, these figures show that the "limit of zero friction" 
discussed previously must be taken literally: It is obtained for such small values of friction 
only that it is practically never realized in nuclear physics. Look for instance at Fig.3.b: 
Already for rj > 0.5 the Dqp is seen to be at least of order T even for T ^ to, whereas in the 
strict limit of zero friction Dgp vanishes identically. 

It has been mentioned that in the computations of regularization of the appropriate 
integral had been achieved by introducing an upper limit in the integral itself, rather than 
to change the integrand applying the Drude regularization. The two methods are compared 
with each other in Fig.4. They are seen to lead to the same dependence on rj but may require 
an optimization of the cut-off parameters, if the comparison is to become more quantitative. 



5 Diffusion coefficients for unstable modes 

The principle which allows one to deduce diffusion coefficients from the fluctuation dissipation 



theorem even for unstable modes has been suggested in 0], see also p9|. In the following we 
like to explore this possibility in more detail with the help of formulas obtained by adapting 
linear response theory to instabilities. All we need to do is to generalize expressions such as 
( ^HD and (|49|) to the case of negative stiffness. 

The basic idea is found in realizing that all the expressions derived and discussed in the 
previous sections are analytic expressions in the transport coefficients, and thus also in the 
stiffness coefficient C. This is true in particular for the general expressions for the diffusion 



coefficients as given by (^Sf ) and (|49|) , together with (p^ . In Fig.5 we present a computation 
of three typical quantities, namely CE^^, Dpp and Dgp, all three normalized to temperature, as 
function of the stiffness parameter C (divided by the inertia M and relative to (7/M)^). The 

18 



result demonstrates the smoothness, or analyticity, across the point C = 0, which actually 
corresponds to the inflection point C(0) = of the static energy. This feature is quite 
interesting, for several reasons. First of all one realizes that the fluctuation in the coordinate 
diverges when C approaches zero from above: S^^ — >• +cxd for C — >• 0"*". Secondly, across this 
point this "fluctuation" gets even negative. Fortunately, for the equations of motion it is not 
this "hypothetical" static fluctuation which matters but the diffusion coefficients. The figure 
nicely demonstrates that all of them remain well behaved at this point implying that the same 
will be true for the dynamical fluctuations! 

On the other hand. Fig. 5 also demonstrates that there are limitations in applying the 
analytical continuation to get meaningful diffusion coefficients. It should be clear that the 
diagonal one ought to be positive. However, depending on the parameters involved, like 
the transport coefficients for average motion or temperature, the Dpp may become zero or 
even negative. This feature can be seen more clearly in Fig.6 where we plot the analytical 
continuations for E^^ and S^^ as function of temperature divided by the bare frequency, i.e. 
of T/zo. The two quantities are normalized such that in the high temperature limit they 
approach the same value, namely 2T /w. Actually, the factor two here is a relict of the 
discussion of the stable case presented in Fig.2, where the same quantities were plotted, and 
where this normalization was convenient at T = 0. Since the curves of Fig.6 just represent the 
analytical continuation of those in Fig.2, it is not surprising that the quantum effects behave 
in the same way as for stable modes when studied as functions of temperature and friction. 

At this stage it may be worth while to refer the reader to the two easy examples discussed 
in the previous chapter, namely the limit of high temperature on one hand and the case of 
zero damping on the other one. Indeed, eq. (|38|) tells one that the relevant diffusion coefficient 
Dpp remains positive for all T, but from eq.(p7D on sees the S^^ to change sign at the point of 
inflection. On the other hand, the T* gets negative below a critical temperature of Tc = Hzu/tt. 
The flgure Fig.6 shows such features to hold true in the general case with the following two 
important modiflcations: a) The Tc is shifted to lower values when rj is increased; b) this shift 
is bigger for the momentum fluctuations. At larger values of t] the latter even show a different 
behavior, in the sense that they become larger than the value given by the high temperature 
limit, which for larger T they approach from above. 

The non-diagonal diffusion coefficient D^p is plotted in Fig. 7 as function of temperature, 
both for positive and negative stiffnesses, and in Fig.S.b by the fully drawn lines as function 
of friction. The values of D^p for positive and negative stiffness differ from each other only 
for small temperatures and small friction, and practically agree with each other already for 
temperatures larger than about O.bhzzi and friction rj > 2. A similar behavior is true for the 
diagonal element, as may be inferred comparing fig. 3. a and Fig.6. At very small T/w the Dgp 
shows an anomalous decrease with T/vu (see Fig. 7 ) or with rj (see Fig.S.b ). This actually 
happens in a regime below T^, the reason of which will be discussed on below. 

Finally, we would like to come back once more to the fluctuation dissipation theorem 
in its general form. We recall that it was this theorem which allowed us to determine the 
diffusion coefficients in terms of the response function of average motion. Furthermore, these 
are exactly the expressions, taken as analytical functions of the stiffness C, which we were 
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able to generalize to unstable modes by analytical continuation from a positive stiffness to a 
negative one. The construction of the special contour C now allows us to perform these steps 
in a more direct fashion, provided one more constraint on the contour is considered. That 
is to say that at instabilities we may express the "equilibrium fluctuations" by the following 
integral: 
@bcoreqinv 

SS = s/,^ooth(i^),CM (50) 

In this integrand there appears another factor having poles in the complex plane, the coth 
which diverges at the Matsubara frequencies u^ = ±n2'7iT/h. As in the stable case, the 
contour C has to cross the imaginary axis in between those which lie closest to the real axis, 
namely ±2TTT/h. This condition actually puts a lower limit on the range of temperatures 
for which such a construction may work, namely T > Tq = /i|u;+|/(27r), with Uj^ being the 
frequency of the unstable mode (see Appendix B, in particular eqs.( |B.9| ) and ( [B.ll| ) and 
the text below them). This follows from the constraint discussed in the context of eq.(^) 
according to which the contour C has to he above the unstable pole. Actually, this Tq coincides 
with the so called "cross over temperature" encountered in dissipative tunneling (see eq.(3.7) 
of (P), which will be commented on in the next chapter. For illustrative purpose we may 
look once more at the limiting case of zero friction for which one gets Tq = hw /{27r) and thus 
2 To = Tc. The feature of the critical temperature Tc being larger than the Tq is generally 
correct. 

Indeed, one may check, with the help of (^l]) , (ji^) and (|43|) , by direct evaluation of (^) 
that this form allows one to obtain expressions hke (^8[) and (^), but now for negative C, 
if only the Drude regularized response functions are used in (EOl) . These features were first 



realized in ^T| and ||T3[ and applied there to treat instabilities in Fermi liquids. 



6 Suramary and Conclusions 

In this paper we have reviewed a method which allows one to describe the dynamics of 
damped motion in collective phase space. The main emphasize was put on the evaluation of 
diffusion coefficients D^,j which govern motion of the second moments both for stable as well 
as for unstable modes. The latter case was handled by a suitable analytic continuation of the 
(quantal) fluctuation dissipation theorem. To this end some elementary properties of linear 
response functions had to be generalized, too. 

From the numerical results it became evident that, when applied to the diffusion coeffi- 
cients, this continuation looses its physical meaning below a critical temperature Tc below 
which even the diagonal element Dpp would become negative. On the other hand, for the 
fluctuation dissipation theorem, as well as for some basic relations of response theory, such 
a continuation is possible even for temperatures down to a Tq < Tc. This Tq was seen to 
be identical to Tq = /i|u;+|/(27r) where u+ = —(0;+)* measures the position of the pole along 
the positive imaginary axis which belongs to the unstable mode. This Tq is thus identical to 
the so called " cross over temperature" [Q one encounters in the phenomenon of " dissipative 
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tunneling" . 

This association to tunneling is by far not accidental. As mentioned above, the equations 
of motion associated to the harmonic case can be exploited for describing global motion within 
a locally harmonic approximation. In the end this means to describe collective dynamics by 
way of a transport equations like the one of Kramers, with the only difference of generalizing 
the classical Einstein relation to the proper quantal form. This implies to have two diffusion 
coefficients Dgp and Dpp instead of the single one Dpp = •yT. Applying this scheme to 
evaluate the rate of decay out of a potential minimum one finds |^ the formula R = JqRk 
known from such model cases for which an application of functional integrals becomes feasible. 
Compared to this latter situation, the derivation sketched here has both an advantage as 
well as a disadvantage. The advantage can be found in the fact that this locally harmonic 
approximation can be applied to cases for which the coupling between the collective degree 
of freedom and the rest is too big to allow for the simplification of linearizing in the "bath 
variables", with the latter simply being represented by oscillators with fixed inertias and 
stiffnesses. On the other hand the description in terms of transport equations of the type 
mentioned is valid only for T > Tc > Tq. In this sense this critical temperature Tc indicates 
a breakdown of the locally harmonic approximation in general. 

In this paper it has been emphasized that the construction of the basic equations of 
motion is possible in self-consistent manor. By this notion we mean that the collective degree 
of freedom is introduced to represent the dynamical microscopic quantity one wants to study. 
This is necessary for all systems which per se exist of identical particles, like it is the case for 
nuclei or the electron gas, to mentioned just two prominent examples. 
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A Extension of linear response theory to unstable 
modes 

In the presence of an unstable mode of frequency a;+ = iT the relations between response 
function have to be modified as compared to the common situation of a stable system. As 
we are going to show in this Appendix, some relations between response function and its 
dissipative part can be retained by making use of deformed paths of integration in the complex 
frequency plane, as encountered in the text for the contour C of (|TBp. 

Causality requires that a physical quantity can be influenced by external fields only from 
their values in the past. The retarded function describing such a response must thus vanish 
for negative times. This may be achieved by an appropriate choice of the contour C^ for that 
integral which defines the Fourier transformation from the frequency representation back to 
time: For negative times, it must be possible to close this integral in the complex frequency 
plane without enclosing any pole of the retarded response function x^(a;). This allows one to 
write: 
@invftresf 

+ 2ze(+t)x"(t) = / ^e— *x^(^) • (A.l) 

Jc+ 2tt 

The contour C+ could for instance be chosen parallel to the real axis if it only lies above the 
a;_|_ = iT (being the pole of x^ having the largest imaginary part). In the case of the advanced 
response, in (|A. 1|) all + signs are to be replaced by — signs, with the contour C_ lying below 
all poles of the advanced function. In the following we will write both versions in combined 
form, using an index X to distinguish between R(etarded) and the A(dvanced). Whenever 
the possible signs ± will appear, the upper one will relate to the retarded, the lower one to 
the advanced response. 

To obtain the equivalent of the well known relation between the frequency representations 
of response function and its dissipative part, one may start with the inversion of ( |A.1| ), 
Oftresf 

X^iuj) = ± I dte'^'2iQ{±t)x"it) \lmuj\ > T , (A.2) 

which is restricted to frequencies whose imaginary part is larger (smaller) than F. In a second 
step, one may combine (|A.2| ) for retarded and advanced response to find an expression for the 
dissipative part of the response function, 
@dissrf 

Under the assumption that retarded and advanced response do not share any common pole 
both contours C± may be deformed to a common one C without changing the value of each 
integral. One only has to ensure that C has all poles of the retarded response on one side and 
those of the advanced one on the other. For a stable system, C can be taken along the real 
axis and the formulas presented here turn to the ones known from the common literature. 
The final step is to insert the dissipative part ( A. 3 ) into ( |A.2| ) and to exchange the order 



of integration for time and frequency. Defining the analytic continuation (with respect to u) 

22 



of the dissipative part of the response function by 
@gendrf 



x"M = ^(x»-x^M) (A.4) 



one obtains the famihar form 
©rftodrf 



;^R(,^) = / ^g!M ioT Imuj > and uj ^ iT (A.5) 

Jc TT U — UJ 

This equation is the key relation for obtaining the integral representation ( ^0]) for the fluctu- 



ation dissipation theorem from the analytic continuation of its representation (^^ . 



B The modified oscillator of the Drude regularization 



In (^51) a frequency dependent friction coefficient was introduced to render some previously 
diverging integrals well behaved. The form chosen there, 'y{uj) = '-f{l—iu!/uj£,)^^, is commonly 
associated to the Drude regularization. Replacing the friction coefficient in the oscillator 
response (|llD by this function the related dissipative part changes to 
©drustrength 

It is apparent that with increasing u the form on right hand side of ( [B.6D drops faster than 
the one with constant friction. Therefore, the Xo^c'^('^) simulates better a correct behavior of 
the true response function at large u, whereas the difference to the mere oscillator response 
function is small in the regions of maximal strength. 

The modified response function has one more pole as given by the new secular equation 
@dispreldru 

C-iuj( + -f) - Mto^ + i—LO^ = (B.7) 

Actually, its solutions depend on two parameters only, as can be seen by introducing x = i-^ 

and rewriting ( p.7|) to 

@dispreldrualt 

-x' + x'- (-^ + —^)x + -^ = (B.8) 

This equation allows for an analytic solution, but it is more convenient to just solve it numer- 
ically. It is, however, instructive to evaluate them in first order in w /wo-, or, less lengthy, to 
solve directly ( B.7 ) perturbatively. Defining rj = 'j/{2Mw) and ( = tu/zoo, one finds for the 



first two solutions 

o;*-^'' ~ ta i -t-i/signC — rj'^ ~ irj + 2r]C 

uj^'^' ^ w { — i/signC — r/2 — ii] + 2riC 
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i^^VsignC -7]'^ - i2r] j 

(B.9) 

--r— #VsignC^^ - i2r] 
1 f/p 



(B.IO) 

It is easily recognized that for r]( —* they turn into the solutions uj± for the common 
oscillator, with cu± = ±JC/M — 'j^/(2M'^) — i'j/(2M). For the third solution one finds: 
©thirdsol 

cu(3) = -iwoil - 2r]() (B.ll) 

A few remarks are in order here. First of all, it is interesting to note that it is not the ratio ( 
itself which matters but the combination with rj. It is intuitively clear that for larger damping 
the Drude frequency must be chosen larger. Secondly, we see that for a large zuc, at least, 
the third pole lies on the negative imaginary axis. 

Let us finally just write down an expression for the full response function which serves as 
starting point for expressing the fluctuations in terms of the \E'— function: 
@chioscdru 

dru ^ Zl UJ + tWp 

where cc;*^*^^ are the solutions of (p.?l ). 

References 

[1] A.O. Caldeira and A.J. Leggett, Ann. Phys. 149 (1983) 374 

[2] P.S. Riseborough, P. Hanggi and U. Weiss, Phys. Rev. A31 (1985) 471 

[3] P. Hanggi, J. Stat. Phys. 42 (1986) 105 

[4] H. Grabert, P. Olschowski and U. Weiss, Phys. Rev.B 36 (1987) 1931 

[5] H. Grabert, P. Schramm and G.-L. Ingold, Phys. Rep. 168 (1988) 115 

[6] H.A. Kramers, Physica 7 (1940) 284 

[7] H. Hofmann, R. Samhammer and G. OckenfuB, Nucl. Phys. A496 (1989) 269 

[8] H. Dekker, Phys.Rev.A 44 (1991) 2314 (for an erratum see: Phys.Rev.E 50 (1996)) 

[9] H. Hofmann and G.-L. Ingold, Phys. Lett. 264B (1991) 253 

[10] D. Bohm and D. Pines, Phys. Rev. 32 (1953) 609 

[11] H. Hofmann and R. Sollacher, Ann. Phys. 184 (1988) 62 

[12] H. Hofmann, "A quantal transport theory for nuclear collective motion: the merits of a 
locally harmonic approximation" , prepared for Physics Reports 

[13] D. Kiderlen and H. Hofmann, Phys. Lett. 332 B (1994) 8 



24 



[14] L.D. Landau und E.M. Lifschitz, Course of Theoretical Physics, Vol. 6, "Fluid Mechanics" 
Pergamon (1963) Oxford 



[15 
[16 

[17 
[18 
[19 

[20 

[21 
[22 
[23 
[24 
[25 
[26 
[27 

[28 

[29 
[30 
[31 



D. Kiderlen, Nucl.Phys. A589 (1995) 320 

H. Hofmann, in: ed. W. von Oertzen, Lecture Notes in Physics 117 Springer- Verl. Hei- 
delberg (1980) 

H. Hofmann, Physica Scripta, T32 (1990) 132 

G.-L. Ingold, Thesis, Univ. Stuttgart 1988 

J. Ankerhold and H. Grabert, Physica A 188 (1992) 568; J. Ankerhold, H. Grabert and 
G.-L. Ingold, Phys. Rev. E 51 (1995) 4267 

P.J. Siemens and A.S. Jensen, "Elements of Nuclei: Many-Body Physics with the Strong 
Interaction", Addison and Wesley, 1987 

H. Heiselberg, C.J. Pethick and D.G. Ravenhall, Ann.Phys. (N.Y.) 223 (1993) 37 

F.A. Ivanyuk, H. Hofmann, V.V. Pashkevich and S. Yamaji, to appear in PRC 

S. Yamaji, F.A. Ivanyuk and H. Hofmann, in Nucl.Phys. A, in press 

D. Kiderlen, H. Hofmann and F.A. Ivanyuk, Nucl.Phys. A550 (1992) 473 

H. Hofmann, F.A. Ivanyuk and S. Yamaji, Nucl. Phys. A 598 (1996) 187 

H. Hofmann, C. Gregoire, R. Lucas and C. Ngo, Z. f. Physik A293 (1979) 229 

M. Abramowitz and LA. Stegun, "Handbook of Mathematical Functions", Dover Pub- 
lic, Inc. New York (1968) 

E.R. Hansen, "A table of series and products", Englewood Cliffs, N.J.: Prentice Hall 
(1975) 

H. Hofmann, G.-L. Ingold and M. Thoma, Phys. Lett.B 317 (1993) 489 

F. Scheuter and H. Hofmann, Nucl. Phys. A394 (1983) 447 

D. Kiderlen, Doctoral Thesis, TU Munich (1994) 



25 



Figure captions 

Fig.l. Equilibrium fluctuations in q (upper panel) and diagonal diffusion coefficient (lower 
panel) as functions of the cut off frequency for various values of damping rate r/ and temper- 
ature T. (The calculations are performed in units with /c^ = 1 and % = 1). 

Fig.2. Equilibrium fluctuations in the coordinate (upper panel) and in the momentum (lower 
panel) as functions of temperature normalized to the frequency w of the oscillator, for different 
values of the damping rate 77. 

Fig.3. Diagonal (upper panel) and off diagonal (lower panel) diffusion coefficient for stable 
(dashed) and unstable (fully drawn) modes, as functions of the damping rate 77, for three 
different temperatures. In the lower panel the curve for T = O.lcu shows O.lDgp/T, as 
indicated by the symbol xO.l. 

Fig.4. Comparison of two regularization schemes: Sharp cut off (thin line) and smooth cut off 
through Drude regularization. The values of the cut off frequencies, u;max and ub respectively, 
were choosen in such a way that in both cases the diffusion coefficient attains about the same 
value. 

Fig.5. Diffusion coefficients (diagonal: thick fully drawn line, off diagonal: dashed line) and 
analytic continuation of the equilibrium fluctuations in q (thin full line) as functions of the 
stiffness C, for fixed values of M and 7. 

Fig.6. Analytic continuations of the (diagonal) equilibrium fluctuations to the case of negative 
stiffness (dashed line for the coordinate, fully drawn lines for the momentum), as functions of 
temperature for different values of friction. The insert shows the momentum fluctuations for 
a larger range of the temperature and the same values of 77. 

Fig. 7. Off diagonal diffusion coefficient as function of temperature for three values of friction, 
for a stable (dashed lines) and an unstable (fully drawn lines) system. The curves ioi rj = 1,10 
show Dqp multiplied by 0.2 and 0.1, respectively. 



26 



1.10 




100 



C0p/G3 



Parameter: r| 

03=0.103 




10 



8 



Parameter: T/G3 

03=0.1030 



>- 



Q 4 







-8 



-10 








Parameter: T/tD 

03=0.1030 



4 6 

11 



<3 







03=1 MeV 




T=0.5 MeV 


^^^^ 


C>0 


^^" " 


-X" 


— Drude- Reg. with 03^=12 MeV 




r* r^O KK> ^1 I'f" /^TT ot /'•'\ OO l\/l/^\/ 


- 


snarp cut on ai oj^^^^— ^u iviev 







0.0 



0.5 



1.0 



1.5 



2.0 



2.5 



11 



! 
i 

Q 



Q 

Q 



iT iT 

W 

o 



2.0 



1.5 



1.0 



0.5 



0.0 



-0.5 




■10 



-8 



-6 





CM/y 



8 



10 



2.0 



1.6 



?„. 1 .2 



o 



B 



0.8 






0.4 












0.2 



0.4 



0.6 



0.8 



T/05 



1.0 



£ 

a 





. 


1 1 1 1 1 1 1 1 1 


1 1 1 1 


0.1 


_ 


0.1 


/ 


^^^,:^^=-^=-^='"^^^''^a2 




1 


— ^ 


0.2 


/ 


^ ""^^ — "^ 






- 


0.3 


/^ 





sign C > 








/ 


Parameter: r| 


- sign C < 










0.4 


- 




03=0.1 03^ 






- 


0.5 


- 




xOA^ 




'■'"^'^lO^ 


- 


n c 


— ' 


1 1 1 


1 1 1 1 1 1 1 1 1 


1 1 1 1 


1 1 


- 



0.0 



0.4 



0.8 1.2 

T/05 



1.6 



2.0 



